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CHAPTER I 


BASIC CONCEPTS 


The concept of inertial navigation is not new. Sir Isaac 
Newton, in his "Mathematical Principles of Natural Philosophy," 
gave the laws of mechanics on which self-contained inertial navi- 
gation is based. At the time of the publication of Newton's 
treatise there were no means available for sensing accelerations 
accurately or for economically determining positions if sensing 
had been possible. It has taken almost three hundred years to 
develop the necessary hardware and software to make self-contained 
inertial navigation a reality. 

A brief description of an inertial navigation system is that 
we have a "black box" which senses both accelerations and mass 
attraction. We eliminate the effect of the mass attraction and 
then use our initial position and velocity, elapsed time, and the 
double integral of the acceleration to determine an updated 
position. In practice this is done many times a second, provid- 
ing a continuous readout of incrementally corrected position. 
Henry described part of the mechanism in the following words: 

A basic requirement of an inertial navigation system 

is that the acceleration vector be known at all times. 

Knowledge of the direction of this vector requires 

either maintaining an acceleration sensitive element 

in a fixed orientation or continuous knowledge of the 


orientation assumed by the acceleration sensitive 
element. 


The system may take any one of the various forms in 
which the acceleration sensitive element or elements 
are maintained relative to the local gravity vector 
(geodetic level) or geocentric level (that is, pointed 
towards the center of mass of the earth). The stable 
element may be fixed with respect to the earth, ro- 
tating around the earth's axis at earth rate, or may 
simply be fixed in inertial space. An azimuth line on 
a level platform may be oriented to point at geodetic 
north or in some other arbitrary direction. (4) 


During the 1950's and 1960's a tremendous effort has been 
devoted to eliminating errors in sensing systems and servo sys- 
tems so that the acceleration sensitive elements may be maintained 
in a known orientation. As stated by Moore: 

The development of inertial navigators has been one of 

the major technological achievements of the twentieth 

century. The outward simplicity of the mechanization 

of these navigation systems and the ease of statement 

of the simple, fundamental scientific principles on 

which they are based mask a major scientific engineer- 

ing and manufacturing development program. It has 

taken the combined efforts of engineering specialists 

in electronics, mechanical, chemical, and electrical 

engineering as well as physicists, mathematicians, 

metallurgists, metrologists, skilled machinists, and 
experienced electronics and electromechanical manufac- 

turing departments working in close coordination to 

develop and produce the apparently simple end product. (6) 

It has become apparent, in searching the literature, that 
some degree of confusion exists about our ability to eliminate the 
effect of mass attraction from the values sensed by the "black 
box."" Statements such as that of Henry, quoted above, equating 
the gravity vector with the normal to geodetic level are common. 
In those cases where some acknowledgment is made that the magni- 
tude and direction of gravity do not agree with simple formulas, 


the effect is considered as random noise in the system. 


It is the purpose of this paper to assess the impact of the 


> & 
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variations of gravity on a vessel, such as a Survey ship, which is 


likely to remain within a limited area for an extended period of 


time, such that the effects of anomalous gravity are systematic 


rather than random, and to suggest some system to eliminate these 


effects from the navigation solution. 


The effects of gravity on the three forms of the system as 


described above by Henry will be investigated to see if the dif- 


ferent forms are affected differently. These three forms, in 


order of increasing complexity, are: 


Ae 


A system in which the orientation of the sensor axes is 
fixed in inertial space, 

A system in which the orientation of the sensor axes 
maintains a fixed relationship to the earth, with the 
acceleration sensitive elements rotating at earth rate 
about an axis which is parallel to the earth's spin axis, 
and 

A system in which the orientation of the sensor axes is 
dependent on the navigator's position, where one axis is 
aligned with the normal to the adopted ellipsoid at that 
position, one axis points north in the plane tangent to 
the ellipsoid, and the third axis completes an orthogonal 


system. 


Three coordinate systems will be employed in this study, all 


of them being right-handed orthogonal systems. The X system is an 


earth-fixed system, with the X* axis parallel to the earth's ro- 


tational axis, positive to the north, the X° axis parallel to the 











Figure l. 


Relationships between coordinate systems 








equatorial plane and the Greenwich meridional plane, positive in 
the direction from the center of the earth to the zero meridian, 
and the X' axis positive in the direction from the center of the 
earth to the ninety degree east meridian. 

The Y system is fixed in inertial space with respect to 
orientation, with Y* = X*, and is related to the X system such 
that, at time t = 0, Y = X% and y' = x’, and at any other time 
there exists an angle of Qt between the respective Y and X axes in 
the sense shown in figure l. 

The Z system is a function of position, and is oriented to 
the ellipsoidal normal such that the Z° axis is parallel to the 
normal and is positive in the direction of the outward normal 
(generally up), the Z axis is positive in the geodetic east di- 
rection, and the Z*? axis is positive toward geodetic north. In 
figure 1, the Z system is shown as being topocentric while the X 
and Y systems are shown as being geocentric. This distinction was 
made in the interest of clarity of the figure only- All three 
systems will be used as geocentric systems so that no translation 
between origins is required. 

Rotation matrices used to rotate the coordinate systems about 
any coordinate axis will be of the form R (aq) where the index "j" 
indicates the axis about which the rotation "q" takes place, with 
the sense of "q'! in accordance with the right-hand rule. Thus 
¥ = Ra (-At)X and X = Rater) Rto)Z where the position vectors are 


defined as columns. 


o 
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The following symbols will be used throughout this study: 


semimajor axis of the ellipsoid 

semiminor axis of the ellipsoid 

eccentricity of the meridian ellipse = (a? - b? )*/a 
gravity vector 

ellipsoidal height 

unit vector along the zero axis of a coordinate system 
unit vector along the one axis of a coordinate system 
unit vector along the two axis of a coordinate system 
Newtonian constant of gravitation 

gravitation vector | 

orthometric height 

radius of curvature in the meridian = a(l-e? )/W* 
radius of curvature in the prime vertical = a/W 

as a subscript, y, indicates navigator 

radius vector 

@ <= sin? »)” 

azimuth 

value of normal gravity 

geoid undulation 

prime vertical component of deflection of the vertical 
geodetic longitude, positive east, measured from XY axis 
astronomic longitude, positive east 

meridional component of deflection of the vertical 
geodetic latitude, positive north 


geocentric latitude, positive north 


g' - astronomic latitude, positive north 


el 
t 


generalized rotation vector 
A  - inertial longitude, positive east, measured from the Y° 
axis =A’ +Qt 
0 - earth angular velocity vector 
QO - first derivative of () with respect to time 
() = second derivative of () with respect to time 
6g - gravity disturbance 
Ag - gravity anomaly 
()° - unit vector in the direction of () 
Vector quantities, other than i, j, k, when written without the 
sign of the vector, ae denote magnitude of the quantity considered. 
There are several quantities used in this study which require 
specific definition in order to avoid a chance of confusion. 
Meridian plane - The plane containing the normal to the 
ellipsoid at a given point and the axis of rotation of the ellips- 
oid. The meridian is the intersection of the meridian plane and 
the ellipsoidal surface for that half of the ellipse which con- 
tains the given point, and is bounded by the intersection of the 
ellipse with the rotational axis. Meridian and meridian plane are 
frequently used interchangeably. 
Prime vertical plane - The plane containing the normal to the 
ellipsoid at a given point which is perpendicular to the meridian 
plane at that point. The prime vertical is the intersection of 


the prime vertical plane and the ellipsoidal surface. It is tan- 


gent to the parallel of latitude at the given point. Prime 





vertical and prime vertical plane are often used interchangeably. 

Latitude, geodetic - » - the angle between the normal to the 
ellipsoid at a given point and a plane perpendicular to the axis 
of rotation of the ellipsoid. 

Latitude, astronomic - »' - The angle between the gravity 
vector at a given point and a plane perpendicular to the mean 
rotational axis of the earth. 

Latitude, navigator's - gy - The value which the inertial 
navigator indicates as being descriptive of one component of its 
position. 

Latitude, geocentric - » = The angle between the radius 
vector to a given point and a plane perpendicular to the axis of 
rotation of the ellipsoid. 

North error - The difference between either the geodetic or 
astronomic latitude and the navigator's latitude, converted to aw 
linear measurement. Whether the error is computed for the geodetic 
or astronomic latitude will be indicated.in the text. 

North error = (gy, - »)M, or North error = (gy - »')M 

Longitude, geodetic - A’. - The angle between the meridian 
plane at a given point and the meridian plane of Greenwich. 

Longitude, astronomic ~ \' - The angle between the astronomic 
meridian plane of a given point, and the astronomic meridian 
plane of Greenwich, where the astronomic meridian plane is defined 
as that plane containing the gravity vector at a given point which 
is parallel to the axis of rotation of the earth. 


Longitude, navigator's - Nw - The value which the inertial 
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navigator indicates as being descriptive of a second component of 
its position, in a direction orthogonal to the navigator's lati- 
tude. 

Bast error - The difference between either the geodetic or 
astronomic longitude and the navigator's longitude, converted to 
a linear measurement. Whether the error is computed for the geo- 
detic or astronomic longitude will be indicated in the text. 

East error = (\y - \)Ncosy, or East error = (Ay = \')Neosy 
Deflection of the vertical - €, n - The angular distance 
between the astronomic position, »y', \', of a point, and the geo- 
detic position, », X\., of the same physical point. The deflection 
is divided into a component in the meridian plane, 4, and a com- 


ponent in the prime vertical plane, j- 


c 


p' = @ 


(A' = \)cosp 


n 

Local vertical rate - The time rate of change of the local 
vertical. This is divided into components in the meridian plane 
and in the prime vertical plane. For an ellipsoid, the meridional 
local vertical rate is the rate of change of the ellipsoidal nor- 
mal in the meridian, or the rate of change of latitude, 0 For 
the geoid, we must consider the latitude rate, plus the rate of 
change of the meridional component of the deflection of the ver- 
tical, or © + E = p'. 

Similarly, the prime vertical local vertical rate for the 


ellipsoid is \cosy, while for the geoid, the rate of change of the 
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deflection of the vertical component must be considered, and the 
prime vertical local vertical rate is Leos ~ n = L' cosy. When 
this term is used in the text it will always refer to the geoidal 
application rather than the ellipsoidal. The meridional local 
vertical rate is also referred to as the astronomic latitude rate, 
while the prime vertical local vertical rate divided by the cosine 
of the latitude is the astronomic longitude rate. 

Gravitation vector - G = The force exerted on a body of unit 


mass because of mass attraction. For a central force field, 


-G.=- SR, where mis the mass of the attracting body. 


Gravity vector -% - The combined force of gravitation and 
the centrifugal force caused by the rotation of the earth, when 
applied to a body of unit mass. 
ge = 6 - Ax (xR) 

Height, ellipsoidal - h - The linear distance, measured along 
a normal to the ellipsoid, from the ellipsoid to a given point. 
Height is positive along the outward: normal. 

Height, orthometric - H - The linear distance, measured along 
the local vertical, from the geoid to a given point. Height is 
positive along the outward normals to the geopotential surfaces. 

Height, navigator’s - h, - The value which the inertial navi- 
gator indicates as being descriptive of the third component of its 
position, in a direction orthogonal to both the navigator's lati- 
tude and the navigator's longitude. 


Height error - The difference between the height and the 








ie! 


navigator's height. The two heights will in general be figured at 
two different horizontal locations. 


Height error = h, - h 








CHAPTER II 
STATIONARY OPERATION OF PERFECT INERTIAL NAVIGATION 


DEVICES IN A NON-NORMAL FIELD 


Conceptually, the simplest form of inertial navigation device 
is one in which the reference frame maintains a fixed orientation 
in inertial space, and which contains accelerometers mounted 
along three mutually perpendicular axes. Given an initial posi- 
tion in space, and an initial velocity, the accelerometer outputs, 
when properly compensated to remove the effects of mass attraction, 
are doubly integrated to provide, at any time, the current posi- 
tion. Although the orientation of these mutually perpendicular 
axes is arbitrary, for ease of consideration we will assume that 
they correspond to the axes of the Y coordinate system as de- 
scribed in Chapter I. 

Equations are developed in Appendix I to permit the determi- 
nation of acceleration in inertial space from the values sensed by 
the errorless accelerometers adopted for this study, under the 
assumption that gravitational forces, which are sensed with accel- 
erations, are normal. This is the assumption used in currently 
operational navigation devices. A similar development, using 
hypothetical accelerations and gravitational field, provides us 
hypothetical values which we consider to be sensed by our error- 


less accelerometers. After completing the general development, 


le 





13 
we consider the specific case of a navigator who is stationary 
relative to the earth at a point where the geoid and an adopted 
reference ellipsoid coincide. The International Ellipsoid was 
chosen for convenience, even though other ellipsoids are presum- 
ably closer to the true size and shape of the geoid. Assumed 
values of gravity anomaly and deflection of the vertical compo- 
nents were assigned to the hypothetical point to which the navi- 
gator was rigidly attached, and these, when inserted into the 
last: equation of Appendix I, provided the apparent motion of the 
navigator under the influence of non-normal gravitation. 

Figures 2, 3, and 4 show the north geodetic error, the east 
geodetic error and the height error caused by components of non- 
normal gravitation. The term geodetic used in connection with the 
word error indicates that the navigator's position is equal to the 
geodetic position plus the geodetic error. Computations were ter- 
minated. when the height error reached + 20 kilometers, which is 
well beyond the limits of certain assumptions made in the mathe- 
matical approximations employed. It should be noted that the ordi- 
nate scale for height error (figure 4) is one-tenth that for the 
north error and east error, and that once the height error started. 
increasing, it was not a slowly changing function as had been 
assumed in obtaining the average values of the gravitation funce 
tion. The results do indicate, however, that this approach is not 
satisfactory. 

If we consider the case where we supplement the inertial nav- 


igation information with an independent determination of the 
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elevation, we are able to limit the rapid propagation of height 
error, although we may have a fairly constant error of the uncer- 
tainty in the geoidal undulation. By using the independent deter- 
mination of height, we not only reduce the height error to zero, 
but also reduce the component of velocity normal to the ellipsoid 
to an independently measured value at each computed position. In 
our specific case, the independently measured normal velocity is 
zero Since the navigator is fixed to the earth. The results of 
using independent height measurements, depicted as before in the 
form of north geodetic error and east geodetic error, are shown in 
figures 5 and 6. 

From a casual inspection of figures 5 and 6 it would appear 
that there is a slight damping of the north error where the deflec- 
tion of the vertical is in the meridian (i-e., & #4 0, 7 = 0), with 
a time dependent increase in the east error, whereas there is a. 
slight damping with time of the east error when the deflection of 
the vertical is in the prime vertical (i.e., &€ = 0, n 4 0), with 
an increase with time in the north error. Figure 7 is a plot of 
the same data.as in figures 5 and 6, but with north geodetic error 
versus east geodetic error, showing that the navigator'’s position 
is oscillating with the Schuler period, and at the same time ro- 
tating as a Foucault pendulum relative to the earth. The example 
was computed for a latitude of forty degrees, and the position of 
the navigator marked at every fifth minute with a small square. 
Although the components of the navigation error are time dependent 


the magnitude of the maximum positional error with respect to the 
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astronomic position remains constant. The magnitude of the maxi- 
mum positional error with respect to the geodetic position varies 
between twice the deflection of the vertical and the square root 
of two times the deflection of the vertical. 

The next step in conceptual complexity of the inertial navi- 
gator, over having the orientation of the reference frame fixed in 
inertial space, is to have the reference frame rotate at the same 
rate as the earth, about an axis parallel to the earth's rotation 
axise This involves only a constant rotational rate about one 
axis, and is independent of the navigator's position on the earth. 
The output of the accelerometers is similar to that discussed 
above, except that the effect of the centripetal acceleration and 
mass attraction at a point maintains a fixed vector relationship 
to the accelerometers, rather than just maintaining a fixed magni- 
tude, and the effects of Coriolis acceleration must be accounted 
for. 

Equations are developed in Appendix II for the case of the 
rotating reference frame, analogous to those developed in Appendix 
I for the inertially fixed orientation of the reference frame. 

The results obtained by using equations II.8 are the same as those 
obtained from using equations I.10, within the limits of accuracy 
of the mathematical approximations made in the derivation of the 
equations. 

The next step in complexity of instrumentation of inertial 
navigation systems is to drive the platform holding the acceler- 


ometers in such a way that the gravity vector always remains 
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normal to the platform. In such a system for surface navigation 
we can eliminate one accelerometer, the one sensing the vertical 
component of acceleration, and measure accelerations only in the 
horizon plane. This system requires that at every point where 
our navigator is likely to be we know the detailed shape of the 
geoid so that we may drive the stable platform at the meridional 
local vertical rate about the east-west axis, at (earth rate plus 
the astronomic longitude rate) times the cosine of the astronomic 
Tatitude about the north-south axis, and at (earth rate plus the 
astronomic longitude rate) times the sine of the astronomic lati- 
tude about the vertical axis. This will provide a continuous 
indication of astronomic latitude and longitude; however, since 
the detailed shape of the geoid is known wherever we navigate, 
according to the requirements of the system, the astronomic posi- 
tion can be converted to geodetic and the problem is solved. 

A variety of things have conspired to militate against the 
successful operation of such a system. One of the problems is 
that of sheer complexity in trying to describe the geoid, to com- 
pute radii of curvature, and to drive the stable platform at the 
correct rates. Another point which should not be overlooked is 
our near-total lack of knowledge of the shape of the geoid in the 
ocean areas; however, for the purpose of this investigation we 
shall assume that we do know the slope of the geoid with respect 
to the adopted ellipsoid at our geodetic position. 

The practical implementation of this type system has taken 


several different forms, all of them assuming that gravity is 
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normal and that the earth is an ellipsoid, and that the stable 
platform is "locally level."' That some of these assumptions are 
not fulfilled is obvious, but the degree of lack of fulfillment 
may be illustrated by the description of the inertial navigation 
equipment used by the U.S.S. Skate in her transit of the Arctic 
Ocean as given by Dozier: 

The N6 platform subsystem mechanization to be described 

may be regarded as theoretically exact and error-free 

if it moves on a perfectly spherical surface and there 

are no imperfections in the components, such as gyro- 

scopic drifts. (2) 

Under the assumption of the spherical earth it was possible 
to divorce platform stabilization from computer operation, per- 
mitting the continued operation of the system in the event of a. 
computer malfunction, which is certainly a very desirable feature; 
however, this system does introduce certain errors into the com- 
puted position because of the fact that the platform is not always 
level. Additionally, in other developments, as a first order 
approximation, height has traditionally been considered to lie 
along the direction the platform is being maintained perpendicular 
to, whether it be the geocentric radius vector, the gravitation 
vector, the gravity vector or the normal to the ellipsoid. 

In order to get away from the various approximations made in 
the literature, the kinematic equations of motion for an inertial 
navigator on an ellipsoid of revolution are derived in Appendix 
III, being presented as equations III.8. Following that deriva- 
tion, a treatment analogous to that in Appendixes I and II for the 


case of a stationary navigator in a non-normal field is carried 
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out, resulting in the average accelerations over a time step shown 
in equations III.16 and III.17. 

When these equations were evaluated with the conditions of 
deflection of the vertical and gravity anomalies used in the two 
cases above, and the results doubly integrated, the errors were 
the same as before, indicating that the alignment of the platform 
axes with the normal to the ellipsoid at the navigator's position 
rather than at the geodetic position was a source of error of 
second order, and, under the conditions of non-normality of field 
which may reasonably be expected, may be disregarded. Although 
the value of position geodetic error approaches one mile at the 
worst condition, the difference in error between the system ro- 
tating with the earth and the system seeking the ellipsoidal nor- 
mal is a maximum of 0.51 meters. 

There are many variations possible in the design of the sen- 
sor orientation, but they are only minor changes from the three 
cases considered here. They may reasonably be expected to react 
in a manner similar to those investigated, so a detailed investi- 


gation of each type will not be conducted. 








CHAPTER III 
OPERATION OF A PERFECT INERTIAL NAVIGATION DEVICE 


TRANSITING A NON-NORMAL FORCE FIELD 


In the previous chapter the operation of various basic types 
of inertial navigation devices was investigated for the restric- 
tive case of being stationary with respect to the surface of the 
earth. Since each of the three basic types behaved in the same 
manner, within the limitations of the mathematical methods em- 
ployed, we will now concentrate on the type which is representa- 
tive of the majority of instruments now available, i.e., the case 
wherein one axis attempts to remain aligned with the normal to 
the reference ellipsoid. 

Equations are developed in Appendix IV for the evaluation of 
the acceleration in latitude and longitude, accepting the fact 
that an independent means must be used to determine the height of 
the navigator if this coordinate is desired. The navigation 
problem undertaken using these equations considered the case of a 
surface vessel, which maintained a constant speed of twenty knots 
and a constant heading of thirty nine and one-half degrees. The 
fact that the vessel remained on the surface constrained the 
orthometric height to be zero, and the navigator's estisave of 


the ellipsoidal height rate to be zero. 
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A. geoidal undulation varying from zero to one hundred meters 
was considered, with a complete undulation cycle Occurring in ten 
degrees of latitude, and fourteen and two-tenths degrees of longi- 
tude. Although the height coordinate, measured from the geoid, 
remained zero, the height rate with respect to the ellipsoid was 
not zero, because of the geoidal undulations. Figures 8 through 
22 show the values of the components of the deflection of the 
vertical, & and n, as a function of time. Figures 8, 9, and 10 
also show the error of position, relative to the astronomic posi- 
tione In this simulation, which ran over a navigation period of 
39 hours, the maximum error relative to the astronomic position 
was 5.6 meters, occurring just after the twenty-fifth hour. This 
figure can be compared to the error with respect to the geodetic 
position of as much as 1750 meters (0.95 miles). 

The results of this test can be generalized as indicating 
that, while assuming that navigation is taking place on an ellips- 
oid and that the force of gravity is everywhere normal to this 
surface, the inertial navigator (mechanically perfect), when ini- 
tially set to proper position and velocity, will track the actual 
direction of gravity on the geoidal surface to a very high degree 
of accuracy, and will provide the angle between the gravity vector 
and the equatorial plane as latitude and the angle between the 
plane which contains the gravity vector and which is parallel to 
the rotational axis of the earth and the Greenwich meridional 
plane as the longitude. Thus, if a map of the deflection compo- 


nents were available with an argument of astronomic position, it 
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would be a simple matter to transform the position so obtained to 
the geodetic position. The effect of changing the gravity anomaly 
from plus to minus was investigated for a period of 1000 minutes, 
with a maximum difference of position of eleven centimeters which 
occurred at 825 minutes when the value of the anomaly was + 46.3 
milligals. Although the anomaly value went as high as + 76 
milligals during this simulation, the position differences between 
the two conditions decreased after 825 minutes. 

Although this study considers only a perfect inertial system 
and blithely ignores all the problems that have been plaguing 
engineers for years, we know that position errors do creep into 
the inertial navigator over a period of time. For this reason, it 
is necessary to reset the equipment periodically to the correct 
position and velocity. Let us assume that some external naviga- 
tion system, such as the Navy Navigational Satellite System de- 
scribed by Russell (9), will provide errorless positioning infor- 
mation in an absolute system, and additionally will provide error- 
less information on the component of velocity in the north direc- 
tion. We will use this external source of errorless information 
every eight hours for resetting our inertial navigator to geodetic 
position and actual velocity north. The results, figures ll, 12 
and 13, indicate that a positional error of 3-2 kilometers (1). 78 
miles) has been reached, with respect to the geodetic position, 
within 45 minutes of the second reset. In view of the fact that 
the navigator attempts to track the local vertical, let us consid- 


er the error with respect to the astronomic position, but under 
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the same resetting conditions. As shown in figures 14, 15 and 16, 
the maximum error is reduced to 1.66 kilometers (0.89 miles), 
which 1s still considered excessive for any type of survey, or for 
any other purpose requiring precise navigation. 

In order to exploit the fact that the inertial navigator 
seeks the local vertical, a simulated run was made wherein the 
position was reset not to the geodetic position, but to the astro- 
nomic position, while the velocity north, as obtained without er- 
ror from the Navy Navigational Satellite System, was used without 
change. The results are shéwn in figures 17, 18 and 19. The 
reader is cautioned about the change of ordinate scale between the 
different figures. Although figures 17, 18 and 19 look, at first 
glance, as discouraging as figures 11 through 16, the maximum 
error, with respect to the astronomic position, is less than 52 
meters (0.028 miles). This error is of a reasonable magnitude, and 
should be acceptable for most purposes; however, for purposes of 
special surveys there are instances in which this error may 
be excessive. 

In navigating through an area where the deflection of the 
vertical is changing, the astronomic latitude rate of the naviga- 
tor consists of two independent components. One of these is, of 
course, the velocity of the vessel over the surface of the ellips- 
oid, which, when divided by the appropriate radius of curvature, 
becomes the angular rate of the vessel. This item, for the lati- 
tude rate, was accounted for in the analysis above. The second 


item which must be considered is the time rate of change of the 
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deflection components. In the present case, in which we consider 


external data affecting only position and latitude rate, we set 


- Von  d 
fr =m t & = + & 


gy is the astronomic latitude rate of the navigator 
Vn is the component of velocity in the north direction as 
provided by an external data source 
M, is the meridional radius of curvature for the navigator's 
latitude 
& is the time rate of change of the meridional component 
of the deflection of the vertical. 
Figures 20, 21 and 22 show the results of this simulation. 
The maximum error introduced in this case, with respect to the 
astronomic position, was less than 12 meters (0.0065 miles), which 
is considered completely acceptable for even the most stringent 
navigation requirements. It should be noted that the geoidal un- 
dulation was not considered in determining the astronomic latitude 
rate of the navigator. Presumably better results could have been 
obtained if this had been accounted for, but since the magnitude of 
the undulation is so small compared to the ellipsoidal radius of 
curvature, which ranges between 6.3 and 6.4 megameters, it was 
considered an unnecessary complication to include it here. 
The environment in this test case is considered to be quite 
gevere, in that an undulation of 100 meters is periodic in less 


than 850 miles. For a specific case of navigation, considering 
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the slope of the geoid in the area of interest and the maximum 
speed of the vessel to be navigated, it may well be that the in- 
clusion of E is not warranted and that the velocity north divided. 
by the meridional radius of curvature is a sufficiently accurate 
approximation of the astronomic latitude rate to provide the navi- 
gational accuracy required. In any event, it is considered that 
if E is to be considered, it should be determined for the geodetic 
position from a plot of deflection components, and entered: man- 
ually into the navigational computer rather than utilizing storage 
in the computer to describe the geoid and to compute the deflec- 
tion rate within the computer. 

We have considered only the use of the Navy Navigational 
Satellite System to provide external data for the inertial navi- 
gator. Other systems, such as bottom navigation, could provide 
not eai'y the information considered above, but also the velocity 
in the east direction, which would be handled in a manner analo- 
gous to the velocity north to further reduce the error after reset 
of the inertial navigator. 

In order to navigate and obtain the geodetic position at all 
times, it would be necessary to compute the components of gravity 
in the horizontal plane at each step of the navigational position 
updating which takes place several times a second in most naviga- 
tors. This would require extensive computer time and increase the 
storage requirements greatly.. By the simple expedient of reset- 
ting the navigator's position to the astronomic rather than to the 


geodetic position, and by inserting the astronomic latitude and. 
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longitude rates (which include the rates of change of the deflec- 
tion components), we are able to accomplish precise navigation 
without introducing an increase in the time required per itera- 
tion, or increasing the size of the navigational computer. In 
order to use the information thus produced: by the navigator it is 
necessary to consider the difference between the local vertical 
and the ellipsoidal normal, but this can be done by graphical 
means or in another computer, and done only when required, rather 
than at each step of the iterative navigational problem. 

We have considered only the case where the externally sup- 
plied information is errorless. Considering the pendulum type 
effect of positional errors, we can conclude that an error in the 
position to which the inertial navigator is reset, regardless of 
cause, will cause the half-amplitude of the error shown in the 
figures (17 through 22) to be increased’ by the amount of the 
error. Since the Foucault rotation must be considered, the half- 
amplitude of the individual components cannot be increased accord- 
ing to the north error or east error, as appropriate, but the mag- 
nitude of the horizontal error must be added to the magnitude of 
the indicated error. The positional error at reset can be caused 
by either incorrect data from the external navigation device, or 
by errors in the determination of the deflection of the vertical. 
The error which is indicated in figures 17 through 22 is the 
result of velocity errors at the time of reset, primarily, and the 
result of the difference between the mathematical model we have 


used for the earth and the true earth shape secondarily. If we 


4? 


consider the maximum error shown in figure 21 (12 meters) to be 
caused entirely by a velocity error at reset, it represents 0.03 
knots error. An error of 0.1 knots (0.0515 meters per second) 


would result in a positional error of 41.5 meters half-amplitude, 


with a Schuler period. 
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CHAPTER IV 
MAGNITUDE OF DEFLECTION VALUES WHICH MIGHT 


REASONABLY BE ENCOUNTERED AT SEA 


Now that we have seen the effect of the deflection of the 
vertical on the operation of a perfect inertial navigator, we may 
logically turn to the question of the magnitude of deflection of 
the: vertical values which might reasonably be encountered at sea. 
If a complete gravity survey of the world were available it would 
be possible to map the deflections for any area of interest to see 
whether or not they were large enough to be concerned about, and 
to make compensation for when resetting the inertial navigator. 
In the absence of an adequate world-wide survey, it is necessary 
to turn to other means to determine whether or not we need an ac- 
curate determination of the deflection of the vertical in the 
ocean arease 

Rice (8), in his study of the area 300 kilometers by 1800 
kilometers reaching from Cape Kennedy, Florida, to the southwest 
corner of Puerto Rico, which is largely an ocean area, obtained 
deflections as great as 56.7". This area can be considered as 
extreme, since on one side of it there are extensive banks which 
are only a few fathoms deep, with the other side plunging rapidly 


to several thousand fathoms, and we would be inclined to avoid it 
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if we wanted to use inertial navigation without compensating for 
the deflection of the vertical. Needless to say, such a value is 
not typical of the entire area. 

Let us consider bathymetric features which are covered by a 
minimum of 1000 fathoms of water, so that the first 1000 fathoms 
below the navigating unit is homogeneous. There are features in 
the ocean areas which are of such an extent that they may be con- 
sidered to be of infinite length, and may be approximated as being 
made of a combination of several bodies which are rectangular in 
cross section. Heiland (3) has given formulas to compute the 


attraction of infinite rectangles as follows: 


Ye ¥ 
ter = eae E] el 


a) 5° 


- tan 3| = a{tan™ 2 - tan 


where o is the density difference between the material of the dis- 


Olt 


+ D {tan 


turbing mass and that of the homogeneous material surrounding it. 
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The above equation may be rewritten for the horizontal rather 
than the vertical component of the change in gravity created by 
the disturbing mass, in that horizontal direction perpendicular to 


the axis of the disturbing mass, as 


t 
- sang] - alten B= tan 
tan Hq h/tan tan 7% 


We will assume that we-are working on a plane, so that Ag, is 


perpendicular to Bg, making it possible to represent the deflection 


by a Three cases have been investigated which are considered: 
as typical of the types of bathymetric features which may be found 
in mid ocean, and near which operations may presumably be con- 


ducted. \ 


CASE I. The Mid-Atlantic Ridge has been approximated by 


1000 fathoms 







70 miles 


1500 fathoms 


100 miles 


140 miles 


Figure 23. Approximation of Mid-Atlantic Ridge 


™ 2000 fathoms 


2500 fathoms 


three blocks of 500 fathom thickness, with widths of 70, 100 and. 
130 miles, spanning the depths between 1000 and 2500 fathoms, as 


shown in figure 23. 
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CASE II. The ridge southeast of the Japanese Islands, near 
» = 32°, &\ = 139°, is similar to the Mid-Atlantic Ridge, but rises 


out of a more shallow ocean. It has been approximated by two 


1000 fathoms 






77.5 miles 


112.5 miles 


Figure 24. Approximation of ridge south of Japan 


1500 fathoms 


2000 fathoms 


blocks of 500 fathom thickness, with widths of 77-5 and 112.5 
miles, spanning the depths between 1000 and 2000 fathoms, as 
shown in figure 24. 

CASE III. The northern part of the Marianas Trench, near 
» = 32°, % = 142°, drops out of an ocean 3500 fathoms deep to a 


depth of more than 5000 fathoms. It has been approximated by 


fathoms 








4000 fathoms 
fathoms 


5000 fathoms 


Figure 25. Approximation of Marianas Trench 


three blocks, each 500 fathoms thick, with widths of 43, 27 and 
13 miles, spanning the depths from 4500 to 5000 fathoms, as shown 
in figure 25. It is by far the smallest feature of the three. 
considered, and is covered by three and one-half times as thick a 


layer of homogeneous matter; therefore, we might logically expect 








we 


a much smaller deflection in Case III than in Cases I and II. 

The magnitude of the deflection in each case is shown in 
figures 26, 27 and 28, as a function of the distance from the 
centerline of the disturbing masses. The maximum deflection is, 
of course, from Case I, being 52.9" at 50 miles from the center- 
line, with a deflection of over 11" remaining 2° therefrom. In 
Case II a maximum deflection of 24.7" was found 56 miles from the 
centerline, with over 7" remaining 2° from the disturbance. Case 
III, with the minimum effect, had a maximum deflection of 18.7" 
at 20 miles from the centerline, diminishing to 2.9" at a distance 
oo. 

Ordinarily, a deflection value of less that 3" of arc would 
be considered negligible; however, if we reset the inertial navi- 
gator to the true position without considering the deflection of 
the vertical (i.e., to the geodetic position), after 42 minutes 
we would have a position geodetic error of 0.1 miles. The im- 
portance of this error would depend on the mission of the vessel 
using the navigator, and, while perfectly acceptable for many 
missions, this error, in conjunction with the error build-up in 
practical navigators, could be the cause of the failure of some 
other mission. 

It should be noted that the bathymetry used in Case II and 
Case III occur with a longitude difference of only oe or a spher- 
ical separation of about Daoee If the navigator were between the 
two disturbances, the effects would be additive in magnitude. In 


each of these cases we have considered the effect of only one 
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disturbing feature, and have completely ignored the compensating 
effect of any isostatic compensation, which would be of increasing 
importance as the distance from the disturbance increased. 

If we consider the world as a whole, rather than isolated 
features, we may work with the gravitational potential as deter- 
mined from satellites, and see what the undulations of an equi- 
potential surface derived from satellite observations are with 
respect to an ellipsoid of rotation. First, let us recognize that 
the potential determined from satellites at altitudes of one or 
more megameters cannot adequately describe the potential on the 
surface of the earth, so the undulations determined are not those 
of the geoid. They represent, rather, the undulations of an equi- 
potential surface of a harmonic potential function in which the 
coefficients of the higher order terms have been inadequately 
determined because of the attenuating effect of the distance be- 
tween the masses creating the potential and the altitude at which 
the effects are measured, in conjunction with the inadequacy of 
the tracking equipment to detect small movements. Additionally, 
the insufficiency of numbers of tracking equipments to cover the 
entire trajectories of these satellites gives a fictitious poten- 
tial function. 

In spite of the shortcomings of the potential function which 
we obtain from satellite observations when we want to apply that 
function at the surface of the earth, the undulations obtained 
from the satellite derived potential should correspond in a gross 


way to those obtained from the gravimetric method, and the slope 
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of the undulating surface with respect to an approximating ellips- 
oid should agree with the average deflections of the vertical ob- 
tained by the gravimetric method when compared in an area of great 
enough extent, as determined by the degree to which the satellite 
potential is accurately computed. Thus, if we were to compute the 
slope of the satellite potential surface relative to the ellipsoid, 
we would not expect to agree in any measure with the deflection of 
the vertical as obtained from gravimetry; however, if we obtain 
the average slope of this surface over an extended distance, we 
would expect fair agreement with the average of the deflections of 
the vertical over the same distance. 

The gravitational field representation developed by Anderle 
(1), which is developed through the seventh degree, sixth order, 
was converted to undulations as outlined by Nicolaides and Macom- 


ber (7). The formula for the undulations is 


7 n 
C = |) 3 P ™(sing)|C, cosm. + s_,oinm,| 


n=2 m=O 


- JeP3 (sing) - 14 (0ind)| 


where 
9 is the radius vector 


Ps are Legendre Polynomials 


Pp” are Associate Legendre Functions 


Cam and oom are coefficients in the gravitational potential 


expansion, and 








58 


Jg and Jg apply to the reference ellipsoid 


The Ja and J4 used were for a reference ellipsoid flattening 
of 1/298.3. From these undulations, deflections were computed 


from the equations 


= AG 
Sy A 


= ee 
1 = WN cosp Av 


The results of this analysis showed that in the North Atlantic 
Ocean, in a region bounded on the south by » = 30°, between the 
meridians of - 30° and - 4? , and on the north by » = 70°, between 
the meridians of - 30° and - 80°, the average slope is on the order 
of seven and one-half seconds of arc. This is, admittedly, the 
largest ocean area which shows a systematic slope of this order of 
Magnitude. If a perfect inertial navigator were reset to the geo- 
detic position at a spot where the average deflection were appli- 
cable, the position geodetic error after a period of 42 minutes 
would be one-quarter mile. Since the average error is of the 
order of seven and one-half seconds, it is reasonable to expect 
some areas where the deflection is larger, resulting in position 
geodetic errors larger than one-fourth mile if the navigator were: 
reset without compensating for the effect of the deflection of the 
local vertical. 

There is nothing in this chapter which should be interpreted. 
to indicate that deflections of the vertical of any given magni- 


tude do exist in the ocean areas. The purpose of this chapter is 
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to show that one cannot reasonably exclude the possibility that 

deflections exist in deep ocean areas, which are of such a magni- 
tude that they could, if ignored, introduce sufficient error into 
an inertial navigator to cause the failure of some mission which 


requires precise navigation for its successful accomplishment. 








CHAPTER V 
CONCLUSIONS AND RECOMMENDATIONS 

A. mechanically perfect inertial navigator is extremely sensi- 
tive to the effects of the magnitude of the gravity vector if an 
attempt is made to navigate along that vector. A slight error in 
the magnitude of gravity will cause an unbounded error to develop 
in the vertical direction. This unbounded error source makes the 
inertial navigator totally unsuited for the navigation task in 
this direction. If we consider a case in which there is an in- 
dependent input for height information, and the inertial navigator 
is used only for navigation in the horizontal plane, the magnitude 
of gravity is of no practical importance. 

When we use a mathematical model for our navigation problem 
which is based! on an ellipsoid of rotation for the figure of the 
earth, with position defined by the normal to that ellipsoid, we 
find that the inertial navigator seeks the local vertical rather 
than the normal to the fictitious reference ellipsoid. If we set 
the navigator to the geodetic position, the navigator’s position 
will oscillate through twice the angle of the deflection of the 
vertical, with a period of a Schuler pendulum, and the navigator's 
position will lie near a plane rotating about the astronomic 
position with the angular velocity of a Foucault pendulum relative 


to the earth. 
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If we had a perfect inertial navigator, and set the position 
initially to the astronomic position while stationary with respect 
to the earth, then the navigator would provide us continually with 
positions within a few meters of the astronomic positions as we 
traversed an area with varying deflections of the vertical. 

Unfortunately, a mechanically perfect inertial navigator has 
not been developed, so that errors do creep into the system, mak- 
ing it necessary to reset the navigator periodically. Resetting 
to the geodetic position is unsatisfactory since the navigator's 
position will then oscillate with twice the amplitude of the de- 
flection of the vertical. In order to avoid this oscillation and 
at the same time retain the current level of sophistication in the 
computer programs, it is mandatory that the position for resetting 
is that defined as the astronomic position, and that the angular 
velocities we accept at the time of resetting be the astronomic 
latitude rate and the astronomic longitude rate. 

If accurate navigation is required through the use of an 
inertial navigator, it is necessary that some means be available 
for resetting position to that of the astronomic position and 
angular velocity to the astronomic rates, or, alternatively, to 
know the gravity vector precisely, and to compensate for the 
geodetically horizontal components if the navigator is reset to 
the geodetic position and the angular velocity is reset to the 
geodetic latitude and longitude rates. The alternative would 
necessitate a great increase in the computation load for deter- 


mining positions. 
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The ideal case for resetting the navigator, from the stand- 
point of simplicity, is to do so in an area where both the deflec- 
tion of the vertical and the time rate of change of the local 
vertical are zero. It is realized that this condition may be 
hard to achieve operationally. 

It is recommended that, where a valid requirement exists for 
precise navigation, and where an inertial navigator is to satisfy 
that requirement, that a complete gravimetric survey of the area 
be made, and that charts showing the deflection of the vertical be 
prepared, with entering arguments of geodetic position in one case, 
and of astronomic position in the other. This will permit taking 
a geodetic position and converting it to astronomic when resetting 
the navigator, and to compute the time rate of change of the com- 
ponents of the deflection of the vertical which are necessary in 
determining the astronomic latitude and longitude rates, and then 
taking an updated navigator's position and converting it to an 
approximate geodetic position when needed. 

The fact that the effects of non-normal gravity were consid- 
ered on a mechanically perfect inertial navigator should not be 
construed to indicate that the writer feels that the mechanical 
imperfections are unimportant and can be ignored. Rather, he 
feels incompetent to consider them without extensive collabora- 
tion with engineers specializing in the field. Now that the 
effects of non-normal gravity have been investigated while being 
isolated, it would be most appropriate to undertake a study com- 


bining the effects of non-normal gravity and mechanical imper= 
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fections, so that it may be seen where the influence of one over- 
shadows the influence of the other. 

Although the underlying assumption throughout this study has 
been that information on the deflection of the vertical would be 
obtained through gravimetric studies, one should not rule out the 
possibility of developing a sensor which would be able to distin- 
guish between the mass attraction effects, the centrifugal force, 
and the reaction to acceleration. Such a sensor would change 
these conclusions and recommendations completely. Similarly, if 
a sensor could be developed which would continuously analyze the 
gradients of gravity to carry either absolute or relative deflec- 
tions with the navigating vessel, the recommendations would have 


to be modified accordingly. 








APPENDIX I 
DEVELOPMENT OF FORMULAS FOR SENSORS 


FIXED IN INERTIAL SPACE 


The Y coordinate system will be used in this case. The 
accelerometers, which sense along the coordinate axes, respond 
to mass attraction and accelerations. Streeter (10) has shown 
that iPetleratione due to rotation about the sun can be neglected 
and that the mass attractions of the matter outside the earth are 
negligible. Letting As represent the value sensed by the accel- 


erometers, we may write 

re = Ay, =e Tel 
= Ane Be -@ 

where Ay, = the acceleration of the navigator with respect to 


inertial space, 


tne = the acceleration of the navigator with respect to 
point B, said point being fixed to the earth on the 
surface of the reference ellipsoid, and 

is = the acceleration of point B with respect to inertial 


space. 


Since point B is fixed to the reference ellipsoid, the 


acceleration is merely centripetal: 
Ze, = Ox(AxR) = Ax(Axn, ) 
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where Ng is the vector from the spin axis to point Bwhichis nor- 


mal to the reference ellipsoid at point B. 
Ase= Ox(OxN,) +0 Raguae G i 
In the problem of practical navigation, the only quantity 
measured is the sensed value provided by the accelerometer. 
Since the value for the gravitational vector is imperfectly 
known, ft is approximated by "normal" gravitation, providing, 


in general, a false indication of the acceleration of the navi- 


gator, Rn» - 


Rye = ie ~ G + Ge I.3 


Normal gravitation has been developed (Hirvonen and Moritz, 
(5)) in terms of the geocentric radius, R, the geocentric lati- 
tude, o, and the longitude, \., once the defining parameters of 
the gravitational reference ellipsoid are known. For the Inter- 
national Ellipsoid with the International normal formula of 
gravity, when R is expressed in kilometers and G in gals, we 


obtain the following: 


G, = - 1000 F + 0.835 888 FP = 2.507 664 Fcos29 + 
+ 0.005 118 F + 0.003 937 Fecos2g - 
- 0.013 778 F’cos* 29 - 
~ 0.000 007 F* + 0.000 054 F*cos2y + 


+ 0.000 023 F*cos* 29 - 0.000 085 F*cos® 29, 


where F = So pent and G, is the component of normal gravi- 
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tation along the geocentric radius vector. The component in the 


meridional plane orthogonal to the radius vector, Ge. is 


Gs = sing cos | - 3.343 551 FP 4 
+ 0-003 149 F - 0.022 045 Fcos2y + 
+ 0.000 031 F* + 0.000 026 F*cos29 - 


- 0.000 145 F* cos? 26 | 


The components of gravitation along the earth-fixed coordi- 


nate axes may then be obtained from 


Gyo = [Gp COS - S5 sing | COBA 


Gy1 = la, cosy = he sing | Sin. 
Gya = G, sing + Gs cosy 


The value of the gravitation vector actually existing at a 
given point is the sum of the normal gravitation vector and the 
gravity disturbance vector, ae. at that point. The gravity dis- 
turbance vector is computed as shown in Hirvonen and Moritz (5), 
and transformed to the coomsinats axes in the same way as the 
normal gravitation vector. 

After we have added the components of the normal gravitation, 
G, and the components of the gravity disturbance vector, 6, it is 
necessary to transform them from the earth-fixed coordinate system, 
X, to the inertial coordinate system, Y, by considering the angle 
through which the earth has turned, Qt. We then put the trans- 


formed components into equation I.3 to obtain the value which the 
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inertial navigator interprets to be his acceleration in inertial 


space, Rye « 

R -|¥° G z “ , + 

Ne + (Gyo + bo - Gyo )eosAt- (Ga + by Gq )sinnt]i 

+| +: (Gyo + 6 yo - Gyo Jeinnt + (Gya + ba - a, cosnt | 3 + Io4 
|? + (Gyo + § ya - Gye )| k 


In the present case, wherein the vehicle to be navigated is a 
surface vessel and remains close to a known equipotential surface 
rather than navigating freely in three dimensions, we may simplify 
the procedure for computing both the normal and the actual gravi- 
tation by using gravity anomalies, Ag, rather than gravity dis- 
turbances, 6g. The magnitude of the actual gravity on the geo-~ 
potential surface W = Wo is normal gravity, computed on the 
spheropotential surface U = Wo, plus the gravity anomaly. The 
direction in which it acts is defined by the angles » + & and 


NK. + n/cosy. 


= 


G= g ~ Ox (AxB) 


Then 


H-4 


Gy = (yy + OQ (Ny + hy) )cospy cos(Ay + Qt) 


eu.$ 


“(yy + O72 (Ny + hy) )cospy sinQhy + Qt) 


= ‘Yn singy k 
Similarly, we may write the expression for the actual gravi-~ 


tation and account for the non=-normality of field. 
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= 


G = - [Cy + Ag)cos(p + &) cos(\ + n/cosp + Qt) + 
7 (N + h)cosy cos(y. + at) | i 
- [Cy + Ag)cos(p + &) sin(i + n/cosp + Qt) + 


Q*(N + h)cosy sin(r. + at) | 5 


md 


- (y + Ag)sin(o + &) 


Substituting these expressions for the actual and the normal 
gravitation into equation I.3, we obtain an expression for the 


value of the acceleration used in the navigation problem. 
Bue = i + (y + Ag)cos(y + &) cos(y + n/cosp + Qt) + 
+ 2° (N + h)cosp cos(\. + Qt) - 
- (vy + O39 (Ny + hy) )cospy cos(ly + at) | i + 
ais + (y + Ag)cos(y + &) sin(i + n/cosy + Qt) + 
+ 2° (N + h)cosp sin(y + Mt) - Io5 
- (Yq + Q* (Ny + by )) cosy sin(Ay + at)| F = 


+|¥" + (y + Ag)sin(y + & - Yu singy | i | 


The gravity formula used in this study is the formula for 
the International Ellipsoid, modified to include a term accounting 
for free-air reduction where the height is positive, and a term 


accounting for the free-air reduction plus the effect of salt 
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water in those cases where the height is negative. The formulas 


used are; 


aco For h positive - 


Y = 978.049 (1 + 0.0052884 sin*» - 0.0000059 sin*® 29) - 


- 0.3086 h x 107° centimeters/sec* 


be. For h negative - 


Y = 978.049 (1 + 0.0052884 sin? - 0.0000059 sin®2~) - 


- 0.2223 h x 107 centimeters/sec® 


The determination of position in such a rectangular system 


is very straight-forward, and may be developed in the following 


Manner. 


t Vv 


R = K& “ Bt + | f 


0 0 


byl: 


du dv 


R = Ro + Rong t + Rog,t + 


t Vv 


t \y 60 
+ | | Rye du dv + ‘i | Ra du dv 
0 0 O O 


Introducing the components of the initial value of the radius 
vector, the initial velocity of the navigator with respect to the 
ellipsoid and the values from equation I.5, we obtain an expression 


for the radius vector at any time as shown below. 
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t V 


(No + bo )cosgo costo + ant + of Yea du dv 
O 


—— 


t Vv 
R = (No + ho )cosgo sinlo + Yb ya t + i | Yh du dv 
O oO 


t Vv 


(i (l-e? ) + hosing + Y3 ya t + il [ Ye, du dv 
O O 


I.6 
Vv 


t 
+ (GxN, )t + | F Qx (Axl, ) du dv 
O 0 


The last two terms on the right-hand side of the equation 
carry point B in a circle about the earth's spin axis, and can be 
accounted for by a simple rotation about the Y® axis. For an 


incremental time step of At we may write 


cos(Q At) -sin(Q At) O 
Ra(-Q At) = |sin(Q At) cos(Q At) O 
O O 1 


Dropping the subscript yx, as being superfluous once the 
motion of point B with respect to inertial space is accounted for, 
and introducing the symbol Pr Ry ~ Re » we can rewrite equation 
I.6 in terms of incremental time steps, relating the position at 
time t + At to the position at time t as the equations which 


follow indicate. 








yak 


1 a ~ ~ 


or 


y? ¥? Ay? y° 
e = Ra(-M At) |v} + {ay eile] | aw ¢ 
: 
BP) (542) Ya (a) Lay (j) Y) (3) 
” g 9 
Ce eC 


L a |v 1 
te & 
etsy ant 


, ‘ i =p ee aoa oe 
Since we cannot, in practice, measure A err. Ye ¢ 1), we must 


oe 


Bubstitue components of R which we derive from A, and which con- 
tain the effects of non-normal gravitation. Thus we use quanti- 
ties obtained from the values sensed by the accelerometers. In 
this study we shall construct the sensed values by taking true 
acceleration, adding the true (hypothetical) gravitation, and then 
subtracting the normal gravitation. Substituting the quantities 
so derived into equations I.7 we obtain the set of equations shown 


below. 








7e 


y° ve Ay? y° 
y = Rg(-N at) |v + [ay + |v At + 
LY (542) 3) 1AM 3) PPT (5) 


+ (y+ Ag)cos(p + &) cos(\. + n/cosp + Qt) + 
+ zat)? Ya (y + Ag)ecos(p + &) sin(i. + n/cosy + Qt) + 


y? 4+ (y + Ag)sin(o + &) 


+ 2° (N + h)cosp cos(y + At) 


+ (N+ h)cosp sin(y + Qt) + 


- a) 
if 8 


Y” (y + Ag)cos(@ + E) cos(y + 7/cosp + Qt) + 
: Eat)? Y 4 (y + Ag)cos(o + &) sin(iy + n/cosy + Qt) + 


Y®? + (y + Ag)sin(y + &) 


+ 2° (N + h)cosp cos(X. + Qt) 


+ (8 (N + h)cosp sin(X + Qt) = 


(j+1) 
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(yn + CN, hy) comm cos(ky, + Qt) 
ia Eat) (Yn + 9% (Ny + hy )Jcosp, sin(\hy + Qt) zs 


Yn singy 


+(j) 


(Vy + (Ny 


Sa 


hy) )cospy cos(Ay + Qt) 


- E(t)? Gor, +0* (Ny hy) )cospyy sin(A\y + Qt) 


Sa 


Yn Singy G4i+1) 


where At is sufficiently small that linear interpolation between 


the components of gravitation at Susy and en) is permissible. 


In actual navigation, where a computer is allocated to the navi- 
gation task and position is updated several times a second, an 
approach similar to the above is acceptable; however, rather than 
use a small time increment and compute many points per second 
which are not needed in this study, the effect of curvature in 
the quantities to be integrated will be considered. 

For the purposes of this study, the true initial velocity 
and the true acceleration of the navigator with respect to the 
earth will be zero. The latitude and geographic longitude rates 
of the mavigator's position will be the result of variations in the 
gravitation vector, and should be small. Equation I.8 can then 


be written, for use with a larger incremental time step in the 


form shown below. 
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a ¥9 Ay° ¥° 
y? = Ra(-QAt) |v + |Ay? + |v At + 
Wd (441) Mei ¢5) AY) (5) LY (5) 


cos(p + &) Av.Val./cos(i. + n/cosp + at)| 
me 
+ Bw + Ag) | cos( + &) Av.Val./sin(a. + n/cosp + at)| (At)® + 


sin(y + &) 


Av-Val.|cos(. “ at)| 
+ S0P (N + h)cosg Av.Val.|sin(a + at)| (Ate. mae? 


0 J 


Av Val. |(y. + Q7 (Ny + by) )cospy cos(\y + at)| 


- 


= 3 Av.Val. (vn + em CN, + hy ) )cos@n sin(\y + at)| (At)® 


Av-Val. | Ys sings | | 


The most. complex of the average value expressions, (yy + 


MAGN, + by ) Jcosipy (22° | (Ay + Qt) can be written as (f, (y,h) + 


fa (y,h)) fa (op) f& (A + Qt) = Fi (Cy, hf CV. + Qt). Since » and h are 
orthogonal to \. + Qt, we know that the average value of the 
product Fy (y,h)f (AX. + Qt) is the product of the average values. 
The quantities py and hy will change, at least initially, very 


slowly because of the initial conditions imposed. We may take 








(2 


advantage of this fact and write the average value of Fi (y,h) as 


Av.Val. Cn + (P (Ny + by) )eosey = 


rlk 


Cn + ae (Ny + hy )) cocina (3) 


mir 


+ 


er a Of ONG + by ) cosy | (j4+1)? and 


Av.Val. | Yu sings | = B(yusingn ) (5) + BCyw singe) (51) 


For the average value of the longitude term we use the fact 
that the average value of some function, f(x), over an interval 


from a to b equals 


b 
= ‘i f(x) dx. 


a 


Ww +t + AV +Q At 
ili 
An. +Q At 
hoot 


sinx ax = 


= hoe: |- cos(y. + Qt + Ar +Q At) + cos(Xr. + at) = 


= paar [costa + Mt) - cos(r. + At)cos(Ar +MAt) + 


+ sin(y + At)sin(Ar +9 bt) | = 
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1 
= 7O At | cos(y. + Qt)(1 - cos(A\ +QAt)) + 
+ sin(y + Qt)sin(AX\. +2 at)| = 


3 
= cos(y + Mt)/ Se +OAt _ (OA +O AEY Se .|+ 


_ GA +Qbt)® | (AA +O dt) — 


+ eae 3 at)/2 7 = i: 2 | 


Because of the initial conditions imposed, AX~0O, so 


AX +Q At=QAt =z 0.005 for At = 60 seconds. 


(Q At)? = 107? (QAt)*= 107? 
3 4 
(Q bt ~ 107° ate ~ 19722 


Since the error of truncation of a series of this type is 
less than the first neglected term, we are safe in going only to 


the cube of AA. +N At. 


2 
Av.Val. (sin(r + at)| = sin(\. + at) (5/2 — (br. +Q dt) | + 


a 
+ cos(h + At) .) Ad aH: a _ Gn. = At) 


Similarly, 


| 2 
Av.Val. [cos(a + at)| = cos(\. + at). [2 - _ a Het a 


2 
- sin(A + Qt) (5) — ite [2 ; {Oh Ate | 
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Substituting the average values of the longitude and latitude 
functions into equation I.9 we obtain the final equations for 
rectangular coordinates of position in an inertial reference frame 


as shown below. 


re cosQ At -sinNdAt O//¥~ Ay? 
y = |sinQAt cosQAt O/lY, + | AY + 
a 2 | Aya 
“ (Q At)? 
y° cos(p + &) [cos(n. + se + Qt) 2 - rae - 
2 
+ |v At + jcos(y + &) [sin(\. + Suze + Qt) ft - (Q At) + 


= ie sin(y + &) 





Les a ae 
sin(rX + ieee + Qt) 5 





=a 6 OY 
+ cos(ry + sean * Qt) 5 


cos(\ + at) [2 - =e 
+ SOR(N + h)cosp (At)* | sinQr. + Qt) fa -—_— 


¢) 
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Po Q At (Q dt )* 
sin(X + Qt) “Che ja - [eee 


QAt (A At)? 
+ cos(h + at) SBE E RE 7 


[Cvs + 08 (Ny + by ) cosn | (5) . 
= Hae) |[[oye + 07 iy + ay )IcoB—N] (5) + 
Ys eingy | (3) + {Vs sings | (541) 


In +R ( Ngee by ))cosen | ( 541] [cos (\n + at)/2 - 


+ Hn + P(N, + by) )eosen| (5 435] [ein(Ay + at)|2 - 


a 
_ (bry +0 At) |- sin(hy +t) Sd +0 bt [a ; 


a 
~ Wen FORE | costay + at) Mu tObt fy 


_ Gbdry mat At)? ]] 


_ CAAy nat <2 1] 








APPENDIX II 
DEVELOPMENT OF FORMULAS FOR SENSORS 


ROTATING WITH THE HARTH 


In the case where the accelerometers maintain a fixed orien- 
tation with respect to the earth, the X reference system is used, 
where X. is identical to Y at time t = 0. Since the reference 
frame is rotating we must consider the effects of this rotation 
on the measured accelerations. The accelerometers continue to 
sense accelerations in the inertial frame, resolved into compo- 
nents along the X axes, rather than measuring accelerations 


relative to the X reference frame. 
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oo 
=> 


Ty = Ry + Qxp + Ox(fixp) +2fixp, +p, =A 

In our case, where the X and Y systems have a common origin, 
and where the rotation rate of the X system is constant at 2, we 
may write 

A = xix) +2xp, + py tages 

As.in.Appendix I, we may write 


= 


al 


II.2 


Although A is the acceleration in inertial space, it need 
not be sensed along the Y°, ¥° and Y* axes. In the present case 
we will sense along the X°, X* and X*® axes. The gravitation 


vector can be expressed in this system as: 


H-) 


Gy = - (vy + Q° (Ny + by)) cospy cosh 


II .3 


are | 


- (yy + Q° (Ny + by)) cospy sims 


- Yn singy k 
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The expression for actual gravitation will be written in the 
same Manner as Was done in Appendix I, taking advantage of the 
fact that, for this study, the navigational experiment is being 
performed at or near an equipotential surface; therefore, gravity 


anomalies will be used rather than gravity disturbances. 


G = -[(y + bg) cos(p + &) cos(vA + n/cosp) + 
+ 2? Neosy cosh. i 
- [Cy + Ag) cos(» + &) sin(d\ + yn/cosy) + II.4 
+ 2° Neosy sink,| A 


=> 


- {Cy + Ag) sinly + 2) | is 


Putting equations II.1 and II.3 into II.2, we write the 
equation that is used to solve the navigation problem from the 


sensed acceleration: 
As = - OF xo ~ 20xt +h + 
+ (Vy + Q9 (Ny + hy) ) cosy cosh | i 


+|- CF? x +20 + x}, ~ 
II.5 


cr. 


+ (vy + Q° (Ny + by) )cospn sinha | 

+ x2 + Yn singn | k 
To obtain the hypothetical expression for sensed accelera- 
tion, we write the same basic equation, but for the theoretical 
position, velocity and acceleration of the navigated instrument, 
and under the influence of hypothetical gravitation. To obtain 


this result, we put equations II.1 and II.4 into II.2 and obtain the 


results shown below. 
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As = |- (Ao - 20% + xP 4 O* Ncosy cos\. + 


Hl 
a 


+ (y + Ag)cos(» + &) cos(A + n/eosy) | 


fle Q* x? + 2X + 4 N* Ncose Sink. + 


II.6 
+ (y + Ag)cos(» + &) sin(r + n/cose) | 5 + 
+| + (y + Ag)sin(y + 2) | & 
Setting ae = a = O, and equating terms between II.5 and 
II.6, we obtain 

X2 = - 07(X° - %) +20%4 + O*Ncosy cost + 
+ (y + Ag)cos(e~ + &) cos(\. + n/cosy) - II.7a 
- (vy + Q7 (Ny + by) ) cosy coshy 

= - oF (8 - x4) - 20% + 7° Ncosyp sim. + 
+ (y + Ag)cos(o + &) sin(r. + n/cos) - II.7b 
~- (vy + Q° (Ny + by) )cospy sinty 

pee = (y + Ag)sin(y + &) - yusingy II.7e 


In view of the fact that 


¥ = Neosy cosh. X2 = (Ny + hy )cospy cosy 


(Ny + hy)cos@y sina 


X* = Neos sina. pk 


equations II.7 reduce to the following expression. 
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x 


x cos(\. + n/cosy) 


Bi = 20) | eh] + (y + Ag)cos(» + &) | sin(u + n/cosp)| - 


ee 


L 0 | | tan(y - 2) 
Die. 


COB\N 
- Yn COBQn sinkn 
P tangy 


Although the determination of the Xe components involves 
more computations than did the determination of the components of 
7, the integration which follows is still straightforward, and, 
Since the gravitation vector for a point remains fixed with 
respect to the rotating coordinate axes (accelerometer axes), a 
simple mean may be used to obtain the average acoeleration over a 


time increment in place of the more elaborate methods used in 


Appendix I. 


t V 
B= Po + Bot + | | rf du dv 
O 0 


This equation may be written, incrementally, in a manner 


similar to that used in Appendix I, as 


= ~ 1 = ~ 
PN (542) = ONC) * Bu 5)dt + Hat)? | Bu (5) i aw (442) 
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XN XN XN 

x, = | + | xi At + 
3 2 

Nd 441) L XN (5) LN (3) 


| ot 
My 05) + 28 (542) 


iy 2 e ° 
+5 (dt) zs. N (5) - W547) 
L : j 
cos(\. + n/cosy) 
+ S(At)? (y + Ag)cos(o + &) |sin(\ + n/cosp) II.8 
tan(y + &) 
COSA 


- F(Ot)P yw (5 yeo8en (5) sinkny 


tangy (3 ) 


COS\N 
- GCA) yw (5 43)2°80N (447) sinky 


tanyn (j+1) 








APPENDIX II1 
DEVELOPMENT OF FORMULAS FOR SENSORS ORIENTED 


TO THE ELLIPSOIDAL NORMAL 


In the case where the accelerometers maintain an orientation 
with one axis parallel to the ellipsoidal normal, one axis in the 
horizon plane in the meridian, and the third axis in the horizon 
plane in the prime vertical, the Z coordinate system is used pri- 
marily; however, since the accelerations are sensed in inertial 
space, it is necessary to transform repeatedly between the Y 


system and the Z system. The transformation equations used are: 


R, = Ri (-~) RelA) Ry = My Ry 


+cosy cosh +cosp sink +sinyg 
M; = |-sind +cosh 0 


-sing cosh -sing sinh’  +cosy 


As before, As =R-G. Let 


R=B+V 
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cosh 
sind 


O 
Ne* cosy 


cosy cosh 

cosy sink 

sino 

N(l-e*) +h 
PB+ Vv seigic ral 
BB + PB + Vv + Vv 

(OXF )B+ PB + (ax )V + VV 

(Gy xP )B + (GXP)B + (Uy xXxP)B+ PB+ BB+ 

+ yx Vv + Gy xP Vv + Gy xP v + Pv + Pv 

(WX )B + (Wx )V + a x(, XE )B + wy x (wy xV? )V + 
_— . ~— ae — Tee 
- eee VV 2(w_, XB )B + 2(w, xV° )V 

cosy 


0 


-sing 


it 
Oo 
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it 


B Ne* cosp sin2 Ne? } ' 
aw? Y p Sing ~ 


We? (1 ae") sing + 


ve Ne? . . 2(1_.2 
B= - a e* (l-e* )siny sin2y - la cosy 
N & Dis 3 e 2 ea o 
See = e) sing cin2y| yg - BT 2s = =e, sing 
= Ne* ( P hss Ne® 2 “2 
| 35 l-e* )siny sin2y - 1 (l-e* )cosy | g® - 
_ Ne? (1-e? )sing 
We 
} 2 2 , F 
“a AS > ey ted e+ 
” (1.2 4 (123 2(1_.¢2 ‘ 
= et ne?) are sin* 2 + Ee ar ) sin® 29 + ieee as ) 2 cosze | e+ 


Me* (l-e*) oe © 
4 sind + h 
We ~) Y 


*% 


a ea sin* 29 Ne? (1-e* ) cosy 9 é 


Ne? (1-e% ) _— 


+ sine 
ow? : 


6 
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O 
) a 0 
» sink 


p sink + » cosA (Q + i) 


Wy =|-~ cosA + » sink (Q + iA) 
oe 2 -_ 8 @ 
h = eae sin® » og 
M, ( B + VV) = O 


2 - a * 3 ° 
a S e") — me sintyg +1 g? + 


fe be 


a a oe 
i 5 k -sink (Q + i) 
Wax =] 0 O Q+A.] = cosh (Q + i) 


cosh sind 0 I O 





i 75 kK 
Wy xP = 9 sinh -~ cosh Q+X = 
|; cosp cosA cosy sinh sing | 
- siny cosh » = cosy sink (Q + 1) 
= |= sing sinh » + cosy cosA (Q + 
| COBY r 4 
0 
Mp (@p XB) =|Q +X. 
0 
O 
My (By xV° ) =] cosy (Q + \.) 


Y 
2M | (Be XE )B + (oy xv? 7 | = 


0 


a @ @ @ @ 
- 2 ne Gee) sin’ y o(Q + ri) + 2 cosy n(Q + dA) 


2 2 ° °° 
2 ne ne sinp cos” yo + 2 hy 
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III.4 
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i j k - sinh X 
Wp XE — |, 0 x = cosh i 


| cosh sinA 0 |} hye” 


= 


md 


a - 
i. J 


4 


ta 
T 


sink » + cosh 9(2 +i.) -cosA @ + sink 9(Q +X) dX 


| cosp cosh cosp sink sing 


- sing cosh » + sing sind @(Q re.) ee cosy sind X 


=| - sing sind ~ - sing cosA y(Q +X) + cosp cosh xX. 


+ cosy o 


fe) 


My (Wy xV° ) =| = sing o(A s 1) + cose \X., 


oe 








91 


M; | (we XE )B + <a xv? Dv | = 


0 III.5 
= |(N + h) cose X% - (N(1-e?) + h) sing @(A BAD 
|(N(1-e*) + h) g 
i 5 k 
Ws X (Wy XE ) =e 0 Q + re 
= sink (Q +i) cosh (Q +i) O 
- cosA (Q +»)? 
=| ~- sink (Q +) 
O 
- cosy (Q + X)8 
My (wy X(gxB)) =| 0 
sing (Q + x)* 
i 7j k 
Wy X (Wy XV? ) = © sink - @ cosh Q+4kX 
| Sing cosA © - - sing sind © + cosy © 


| - cosp sink (Q + x) + cosy cosA (Q + 8) 
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- cosp cosh mi + Sing sinA e(A + rv.) - 


= = —), @ e e 
Wy X(Wy XV9) =| - cosp sind g® - sing cosh o(M +) - 
fo Sing yo 


- cosy cosA (Q + 1) 


- cosy sink (Q + i)* 


- @ - cos*y (Q + A) 


My (Wy X (Wy xV° )) sing o(Q + 1) 


tt 
f 


+ sing cosy (Q + i) 
M, | ty x (Wy XB) B + Wy X (Dy xV? )v | = 


2SON(@=e?) + h) e= Gi + h)cos* KQ + FoF TEps6 
=| - (N(l-e?) + h) sing 9(Q + i) 


| + (N + h) sing cosy (Q + i? 


Equations III.3 through III.6 provide the components of R 


that go into equation III.2. To determine R as a function of the 
sensed acceleration, it is necessary to obtain the components of 
the gravitation vector also. From considerations similar to those 


of Appendix I, we will derive the gravitation vector from the 
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gravity vector, since we are considering navigation on or near the 


reference ellipsoid in the ideal case. 
~ G+ Ax(x(N + n)V?) = -¥ 


vy + (N + h)cos*y 9 
M, (-G) =| 0 III.7 


L- (N + h) sing cosy 0% | 


o@ ee o"~ a °o ° e 

Ze h (Se + h| ~ (N + h)cos*g 1(20+ 2X) + Y 
—_ oe oe = 3 e e 
A, =| 2) = /(N + h)cosg X= [Base + h| sing 9(Q +A) + 


oe —~ae oe _) a _ °° 
72 [Sa + h | ae uC =) Sx sind g + 2hg + 


+ 2 cosy h(Q + A) 


+ (N + h) sing cosy GIA + du.) 


N(1-e® ) 
2 


In view of the fact that W 


is the radius of curvature 


of the ellipsoid in the meridian, M, we may rewrite the equations 


for the sensed accelerations in the more compact form shown below. 
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Ze h- (M+ hg? - (N + h)cos*o i(2N+ 1) + 
Ze} = 1(N + h)cosy X% - 2(M + h)sing o(Q +i.) + 
73 | (M + h)g Me" 5 °3 oe 
¢ p+ 3 we Sing cosp g + 2hy + 
Lit.é 
+ Y 


+ 2cosm h(Q + i) 


+ (N +h) sing cosy 1(20 + i) 


The above equations, when evaluated at the navigator's 
position, permit the determination of h, X. and * from the 
sensed accelerations, under the assumption of normal gravity. 
These equations are equivalent to the kinematic equations of 
motion presented in standard works on inertial navigation, but 
these equations are in terms of fundamental parameters used by 
geodesists (M, N, e*, », X\., h) rather than in terms which have 
more general applicability, and which are used to fit any situ- 
ation (radius, geocentric latitude, longitude), but which must 
then be transformed to obtain the quantities which the navigator 
desires. 

As in Appendixes I and II, a problem will be considered in 
which a navigation system is used which is initially in perfect 
alignment with the coordinate system at the geodetic position of 
the observer, but in which the gravity vector does not coincide 


with normal gravity. The centrifugal effects given in equation 
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IlIl./7 will remain unchanged, since they are a function of position 


only, but the net gravity effect will change to 


Y + dg 
(ga), = Ri (-~) Re (A)Ra(- A - n/cosy)R (y + &)I0 


10 


| 


Using this relationship, and tacitly assuming that we have no. 
geoidal undulation with which to contend, it is possible to 
rewrite equation III.8 for the sensed accelerations in the hypo- 


thetically true case as 


7 h - (M+ h)g® - (N + h)cos*y 1.(22 + i.) 


@| = |(N + h)cosy % - 2(M + h)sing @(Q me) 


oe ee Me? 7 *o e0e 
| 23 hu + h)p + 3 Var sing cosy gy” + 2hy + 
+ 2cosy h(Q + A) + 


+ (N + h)sing cosy 120 +O 
II1l.9 


sing sin(y + &) + 
+ (y + dg) |e - - - 


| cosp sin(y + &) - 
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+ cosp cos(~ + &) cosh cos(A + n/cosy) + 
- cos(y + &) sind cos(A + n/cosy) + 


- sing cos(y + &) cosA cos(A + n/cosy) - 


+ cosy cos(p + &) sind sin(A + n/cosy) 
+ cos(p + &) cosA sin(A + n/cosy) 


- sing cos(y + &) sind sin(A + n/cosy) 


In the cases which have been considered in the first two 
appendixes, the navigator has been stationary with respect to 
the earth, so that h = © =N=h=qg=ih.= 0. In this situation 


equation III.9 simplifies to 


ze cos(A + 7/cosy) 
Zl = (y + Ag)cos(p + &) My | sin(A + n/cosy) III.10 
pir tan(p + &) | 


If the stable platform were to maintain its orientation with 
respect to the normal to the ellipsoid, and to rotate such that 
the Z* axis remained aligned in the meridian, the above quantity 
would remain constant. Unfortunately, the rate of driving the 
platform about the various axes is a function both of indicated 
velocity and indicated latitude, and the axes themselves are con- 
stantly moving, making the problem difficult. The rotation rate 


about the instantaneous Z axis, which, for the navigator, is 








el 


defined as east regardless of its actual orientation, is Oy e The 
rate about the instantaneous Z* axis is (Q + hn )COBGy 5 and that 
about the instantaneous Z axis is (Q + hn ) 8ingy « Since these 
rotations are not vector quantities, the order in which each 
rotation is accomplished is important; however, these rotations 
are accomplished simultaneously in practice and in infinitesimal 
increments. 

Let us consider that during one time increment, dt, the 


transformation matrix, M,, is modified by the transformation 
(541) - Ro ((Q + Ay ) Singy dt )Rz (=n dt )Rg ((Q + An )COSQN dt )Mr (4) 


This expression would become quite unwieldy when expanded; how- 
ever, if we limit dt such that we may consider the cosines of the 
resultant angles to be one, the sines as the angles in radians, 
and the product of two angles to be negligible, the order of 
multiplication becomes unimportant, and the relationship may be 


approximated by 


Ojai) 7% eM (5) 


where 
+1 +y COSQy dt On dt 

Me = |-Aycosgy, dt +1 Ay singy dt 
=On dt ~Ay singy dt a: 


where Ay =Q + Male At any time t we may write 
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ee 


Be cos(A + n/cosy) 
t 


AA (y + dg) a Mg (v)Mr (4) cos (gp + &) | sin(A + n/cosy) | III.11 
¥= 


<o | 
Zs _ tanly + 16) | 

The evaluation of equation III.11 can be extremely tedious 
because of the small time step required in the evaluation of the 
matrix Me(v). We may rewrite that equation in terms of individual 


rotation matrices in the vollowing manner: 


=e = TT | Ral-oy gat; )Ro( (0 + ln )singy dt, )Ra((Q + 


se iy eosen , at, )] Ri(-)R3g(rX) x Litet2 


x I Ra(-fat ,)| Bx 
where the product Ra(-Adt ;) transforms the gravity vector from 
X space to inertial (Y¥) space. The product Ry (-~)Rg(\) trans- 
forms it to the Z space as defined at time t = 0, and the long 
product which makes up the remainder of the expression transforms 
to Z space as defined at time t, which corresponds to the instan- 
taneous sensor axeSe 


The two rotation matrices 


Ro ((Q + Ny )Singy dt)Re((Q + dy )cosmy dt) 


represent nothing but the component rotations of (Q + ry) about 


the vertical and north axes. These can be replaced by a matrix 
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R2((Q + Xy)dt) placed where the axis of rotation is either X* or 


i. 


(aay 







(Q + hy ) COBO, 
\(Q + Ay) sings 


Thus, equation III.12 may be written 


int [Ra (on gat ) | Ri(-p)Ra(\) x 
i 


pu 
eo 
Hl 


x U | Rat (0 + An ;)dt ,)Ra(-Mat ,)| go 


ll 


TT [ Ra (~ou dt, | Ri (-~) Ra (a) TT [ Rally jdt )] Bx Eases 1.3 


Since the product of rotation matrices about the same axis is 
merely the rotation matrix for the sum of the individual argu- 


ments, equation III.13 may be written as 


t t 


R Ri | ~y at Ry (-~ Ra (i Re f haat Bx 
© 0 


oP 
@ 
ul 


cos(y\. + n/cosy) 


Ri (-oy )Ra(Ay) (y + Ag)cos(p + &) | sin(r + 7/cosy) Tire.4 


tan(y + &) 
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Combining equations III.8 and III.14, we obtain expressions 


for the observed accelerations of the navigator in height, long- 


itude and latitude as follows. 


hy + (M+ h)y o® + 
(N + h)ycosgy tw] = |+ 2(M + h),y Singy oy (2 Pn ee 
(M + h)yon | dF 3Mye*® sin2py g@/2WR - 


+ (N + h)y Ccos* wy Cos Nee) - ‘Yn 


- 2cospy hy (Q + An) + Tae hs 


2hw on - (N + h)ysingy cosgy Ay (22+ Ay) 


cos(\. + n/cosy) 


+ 


Ra (-gy )RaQan ) (CY + Ag)cos(p + &) | sin(r + n/cosy) 


tan(y + &) 


oe 


MW ay = - (y + Abg)cos(» + &) x 


- sint, cos(\. + n/cosp) + cosly sin(\ + yeeee) | f 
Ny + hy )cosgn (5) 


(= sim, cos(\. + n/cosp) + costy sin(y. + uicose? ) 


(Ny + by )cospn i ° 








{Gn bw deinen gui + tw) = cosenbe(Q + dw) 
(Ny + hy ) cosn : 


(3) 


-[Siu+ nae see he) ~ cower SP iGnas aly 
(N, + hy ) cosy 


Collecting terms, and putting hy 


av 
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ye 
Ny COSQn 


, | Musings Oy (Q + Aa) 
Nu COs@n (j+1) 


S = + A zs 
PNay = 2(My + hy) (5) | 


+ A 
i 2(M, + ba) 549) |- 


2 A 5 
7 TH, + ba) | My si n20pn AR | (5) = 


y, « rly + dgdcos(y + &) | /sinl 
— . ~(y +4 gost + sink 
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Sing, 
Singy 


COS Qn 


Singy 
Sing, 


COSYn 


COS\N 
SitAN 


sin( 


COSA, 
sinky 


sin( 


oars 
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(j+1) 


hy 


My singy oy (A + ha) 


Ny COS®y 


cos(@ + &) 


cos(» + &) 
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cos(p + &) 


cos(» + €) 
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—~ 
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= O, we obtain 


PELE GA ILS 
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+ »/cos)- 


+ /cosy)+ 


+ 7/cosy)- 


+ 1/cosy)+ 
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ee Sin2py ge (2IO' + a) | a 
J 
e? . ; o- 
TS aT gap meinen My LMT] 


=| s ‘ sin2wy Lagi ® Gee ka) | 
N N Gey, 


Collecting terms, and putting h, = hy = 0, we have 


ON ay = Bey) |- sing, cos(y + &) cos(An - &\ ~ n/cosp) + 


+ cosy, sinly + 6)] 5) “ 


+ EE | - sing, cos(y + &) cos(Ay 
(j+1) 


\. - n/cosp) + 


+ cosyy sin(y + @)] 543) : rag Y, 


= 2e [ein2y oR WR] C5) - - (ae sin2oy Ay (20 + mw] 5) > 
- 32° |ein2oy oR NR | (j+) -{ ae sin2gy Aw (20 + hw) | Gy 


Once the average accelerations in latitude and longitude are 


computed, the determination of position follows as 


] : Liat)? x 
M( 541) = N45) + A 4) 4t + (ht) Nay 


% 


7 ° i eo 
Tee iy * apo" * BIAty” ay 


The development in the appendix through equation III.15 


would indicate that computation for height would be possible in 
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addition to the computation for latitude and longitude. The 
development was not completed since, as the results of the two 
previous sensor configurations show, a solution for height in an 
imperfectly known gravitational field leads to ridiculous results. 
One of the primary objects of the present orientation of the sen- 
sors is to isolate the latitude and longitude sensors from the 
height sensor, recognizing that it is necessary to supplement the 
inertial data with information from another source for height 


information. 








APPENDIX IV 
DEVELOPMENT OF FORMULAS FOR NAVIGATION ON A. 


HYPOTHETICALLY UNDULATING SURFACE 


In the first three appendixes we have considered the effects 
of anomalous gravity on a navigator which was stationary with 
respect to the surface of the earth. Referring to equation III.9 
we see that by having the navigator stationary we have eliminated 
a large part of the equation, which could, in fact, provide a. 
significant contribution to the sensed value of acceleration. In 
this appendix we shall consider the more realistic case of a: 
vessel moving over the surface of the earth, wherein it is sub- 
ject to a varying effect of gravity, and where all the factors 
which contribute to the sensed value of acceleration must be con- 
sidered. 

Let us consider the case where the geoidal undulations may 


be expressed by 
6 = Go [cos Srth = hen) - 1|[cos ohio =e) - 1] IV.1 


where Ao, AL, ~, and m are the boundaries of the tessera in 
which undulations are to be computed and 4€5 is the value of the 


undulation at the center of the area. Using the relationships 
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we may write analytical expressions for the components of the 


deflection of the vertical in the following form. 


_ bo 2n(y. = ho) _ 2m . 20(@ = go) 
B= [cos CEES mo- we Ce a IV.2 


- 62 Tee 2m . 2u(X. = do) 
1 * Reose [eos fi - ~ ire -o Ce -ast a 


Although § and n are directly dependent on the variations of 
the undulation, €, there is no such rigid dependency between € 
and the gravity anomalies. We may start from the principal equa- 


tion of gravimetric geodesy 


dg = - act - 


VIO 
rales 


and see that since T, the disturbing potential, has not been 
specified in a form with a unique derivative, we are relatively 
free in our choice of the hypothetical distribution of anomalies. 
The tests described in Chapter II seem to indicate that the mag- 
nitude of gravity is relatively unimportant, so long as we have 
an independent means of measSuring height. As a check on the 
importance of dg, we shall let it vary in two different ways, 
then compare the results. In one case we shall assume that Ag 
in milligals equals € in decimeters, while in the other case we 
shall assume the negative of this relationship to hold. 

Let us consider the instant of time at which we are steaming 
at a speed V, and at an azimutha. The relationships which 


relate the rates of change of latitude, longitude and elevation 
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to velocity and azimuth on the surface are 


V_cos(&? + 12)” cosa 
M+ H+ ¢ 


ae V_cos(E* + »* yz sing 
ail (N + H + €)cosy 


V sin(& + n* y2 


ry’o 
i) 


Since, for this problem in navigation, M, N, &, n, 65 9, V; 
H and @w are functions of time, the analytic expressions for Os _ 
and h, which are needed: in equation III.9, get quite lengthy. In 
order to avoid the use of a very small time increment, we will 
establish the simplifying condition that V, H anda are constant; 
however, the fact that H is constant does not imply that h is 
zero, since h must take care of the variations in Go Since the 
establishment of analytic functions for © 5 \. and h is a formidable 
task, even with the simplifying conditions, we will use numerical 
methods to arrive at these quantities. 

We know that 0, d. and h are well-behaved, slowly changing 
functions. Using standard prediction techniques, we will fit a: 
curve to the functions at the four most recent time points, then 
use this curve to predict the value of the function at the time 


point next to be considered, and the point immediately thereafter. 


The difference in the function at | gel) and =) divided by 


the time increment between the two points will give the value, 


to an acceptable order of accuracy, of the time derivative of the 
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function at Las The imposition of an iterative cycle if rigid 
limits on the difference between predicted and computed positions 
and functions are not met insures that we will not use values 
which are significantly in error in any of our equations. 

As a set of equations equivalent to equations I11.10 and 


III.8, we may write 


i cos(y. + n/cosp + Qt) 
Zl = (y + dg)cos(o + &) M |sin(y + n/cosp + Qt) | + 
Z3 | Itan(@ + &) 


h - (M+ H + 6)g" ae Gi. ai > ¢ oat (20 + 2d) 
+ |(N + H + €)cosy X - 2(M + H + 6)sing o(Q + Ome 


| (M + H + Co + 3Me* sing cosy y /W + 2hy + 


+ 2cosy h(Q + X.) = 


+ (N + H + €)sing cosy M20 + ) 


IV.4 


ha ~ (My + hn gh ~ (Nx + hn )cos* on Aes G2 ae) + 
= |(Ny + by )cosgn i - 2(Myn + hn )singn on (Q + >n) + 


(My + by on + 3My e* singn cosen gh/Wn + 2hnon + 
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+ Yn 
+ 2cos~y by(Q + Aw) 


+ (Ny + hy )singy cosgn tll + ee) 


In equation IV.4, hy is measured from the geoid, but the 
Navigator, in applying it, assumes that it is measured from the 
reference ellipsoid and assumes that the geoid coincides with the 
ellipsoid. After going through a derivation similar to that pre- 
sented in Appendix III we arrive at an equation which is the 


equivalent of equation I1I.15. 


ee 


hy’ h- (M+ H+ Ge - 
(Ny + hy )cosgn we = |(N + H + €)coso — 
(My + hy Jon (M+ H + C)o + 3Me*sin29 g@/2Ww + 


(N + H + 6)cos*o Nidameee Le) 


- 2(M + H + €)sing g(Q +X) + 2cosy n(Q+r)| + 


2he + (N + H + €)sing cosy M20 + de) 


a 


IV.5 


cos(k. + n/cosy) 


+ 


Ri (-py )Ra(an ) (Cy + Ag)cos(p + &) |sinly + n/cosp)| + 


tan(y + &) 
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(My + hy )g + (Nn + hy )cos* ox Dei BG + wie) = YN 
+ | 26% + ohydeined o6(Q + Xe) = BeGben bn (OQ +Mee) 


~3Mye®sin2on ga/2Wa - 2hvon - 


- (Ny + hy )singy cos An 62 0 si) 


Expanding this equation we get 


i = ( + Ag cogie + sin(i\. + n/cosy - Aw) + 
- (N + H + eas = 28 fee C )sing g(A + 4) 
Ny + hy )cosgn 
i. 2cosy h(a + d.) + 2(My + hy )singy on (Q - Aa) 
Ny + hy COS ON 
_ 2cosyy hy (Q + dan) 
Ny + Hy )cosgn 


Collecting terms, averaging over a time increment, and 


setting H = hy = hy = O, we obtain the final expression for the 


average longitudinal acceleration. 





7 







7 a Swe ee ee eee | ee —e 


a 
a 7 ns on! ae eke td 6 oo =a 
7 hee eee 


~ 
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oe 1 
Way = cam |e y( 3°" + Ag)cos(» + &) sin(y + 


+ n/cosp - lw) + S(N + €)cos soe 
- (M + €)sing (A + a) + COSY h(Q + i) + 


+ My ain, “oy (AQ + bw] (5) Z 
IvV.6 


a 1 


+ n/cosp - Ay) + =(N + C)cos» AU 
- (M+ C)sing o(Q +A) + cosp h(Q +A) + 


+ M, sin, oy (A ne hw} (542) 


(My + hy on = (¥ + Ag)cos(o + &) cospy | tan(y +e - 
- tangy, cos(\ + n/cos - dn) + 
+. $sin2e| 30° Mg? /W + (N + H + Q)at2n+ 2] + 


+ (M+ H+) + 2h - 2hygy - 


Z Ssin2n | 36? Muon /h + (Ny + by )An(2Q+ hw | 


Collecting terms, averaging over a time increment, and 
setting H = hy = hy = O, we obtain the final expression for the 


average latitudinal acceleration. 








ook 


Nay = EAA (y + Ag)cos(y + &) cosy | tan( re 2) = 


tangy cos(\ + 4/cose - dw) | + 2he © 


+ SEB28/ sme? We + (N+ 6) X(2a+K%)] + 


(M+ GC) - SARE (3m, 0? Gk AR + 


+ 


+ 


taka (20+ ind] 4 
| (5) IV. 


es 
+ ae] ay ( + Ag)cos(p + &) cosgy | tan(e + &) - 


tangy cos(X + 7n/cose - hw )] + 2hy + 


bats [aMe? GW" + Otis ¢) acan + X» | > 


rs 
+ (M+ C)@ - SBE (51 7 on + 
+ 


fie aso in] 
(j+1) 
Once the average accelerations in latitude and longitude are 
computed by means of the above equations, the determination of 
position at any time is the same as that developed in Appendix 


TLE 
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